WIS-94/2/Jan-PH 
hep-ph/9401246 
January 1994 



Phenomenological Constraints on A and A] 



Zoltan Ligeti and Yosef Nir 
Department of Particle Physics 
- - ' Weizniann Institute of Science 

Q>^ . Rehovot 76100, Israel 

a^ 

c . 

c3 ! Abstract 

cn . 

Combining the experimental data on the inclusive decays D -^ X ei^, 
B ^f X ev and B ^f X Tv,we find severe constraints on the A and Ai param- 
^ ■ eters of the Heavy Quark Effective Theory. In particular, we get A < 0.7 GeV. 

t:;;j- . Our constraints further imply rric > 1.43 GeV and mf, > 4.66 GeV. We dis- 

^N I cuss future prospects of improving these bounds and their phenomenological 

Q ' implications. 
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The Heavy Quark Effective Theory (HQET) provides a systematic way of calculating 
l/mg corrections to the model independent results of heavy quark symmetry |l[]. For inclu- 
sive semileptonic B and D decay rates, Chay et al. showed that the heavy quark limit 
(mg — >■ cxd) coincides with the free quark decay model and that there are no corrections 
to this result at order l/mq. The I/itlq corrections can be computed in terms of only two 
additional nonperturbative parameters 0. 

HQET provides much more precise predictions for inclusive semileptonic decays than 
the spectator quark model [Q. One ingredient in the improvement is the incorporation of 
the nonperturbative l/rriQ corrections which are, however, numerically quite small. More 
important is that HQET gives an unambiguous definition of the quark masses , and these 
masses determine the inclusive decay rates. Using the relation between the masses of heavy 
quarks in terms of the HQET parameters A, Ai and A2 (see their definitions below) leads to 
a significant reduction in the theoretical uncertainties 0. 

These HQET parameters are genuinely nonperturbative: at present they can only be 
estimated in models of QCD. In ref. 0] we used such theoretical ranges for the various 
parameters to calculate inclusive semileptonic B decay rates. Here we take the opposite 
approach and investigate what restrictions can be derived from the existing experimental 
data on these parameters. Our constraints help to pin down the physical values of the HQET 
parameters. This will allow testing various model calculations and will lead to significantly 
reduced theoretical errors in the predictions for inclusive heavy hadron decays. 

We now turn to an explicit presentation of the relevant HQET parameters. In the 
framework of HQET, semileptonic decay rates depend on three parameters (besides the 
Fermi constant, CKM angles, and the well-known masses of physical particles): A, Ai, and 
A2. A is given by the matrix element ||^ 

{0\qzv'DrK\M{v)) 

{0\qTK\Miv)) ' ^ ^ 

where M{v) denotes the meson state and h^ denotes the quark field of the effective theory 
with velocity v. This parameter describes, at leading order in l/mg, the mass difference 
between a heavy meson and the heavy quark that it contains, and sets the scale of the l/mq 
expansion. The dimensionful constants Ai and A2 parameterize the matrix elements of 
the kinetic and chromomagnetic operators, respectively, which appear in the Lagrangian of 
HQET at order l/mq: 

Al = ^{M{v)\KiiDfK\M{v)) , (2a) 

A2 = -^{Miv)\^Ka^f,G''''K\M{v)) , (2b) 

where t/p = 3 and dy = —1 for pseudoscalar and vector mesons, respectively. Whereas Ai is 
not renormalized due to reparameterization invariance 0], A2 depends on the scale. 

The dM factor appears because the chromomagnetic operator breaks the heavy quark 
spin symmetry. Consequently, the value of A2 can be easily extracted from the mass splitting 
between the vector and pseudoscalar mesons: 



A2= "'^' ^^^^^ -0.12GeV^ (3) 



m^, — rrif 



We expect this value of A2 to be accurate to within 10%, as a result of the finite b quark 
mass and the experimental uncertainties. There is no similarly simple way to determine Ai 
and A. 

The HQET parameters appear in the expansion of the heavy meson masses in terms of 
the charm and bottom quark masses: 

- Ai + 3A2 / . , 
iTiB = mb + A h . . . , (4a) 

- Ai + 3A2 / .1 \ 
mD=m, + A — + . . . . (4b) 

For each set of values {A, Ai, A2}, eqs.(|) determine rric and nib. The consistency of the heavy 
quark expansion requires that these values of the quark masses are used in the theoretical 
expressions for the decay rates. Then, a precise knowledge of A and Ai would allow accurate 
predictions for heavy hadron decays. 

It was noticed in ref. |^ that semileptonic D decays provide a lower bound on the 
parameter Ai at each particular value of A. At present, there exist three measurements of 
inclusive rates that provide bounds on A and Ai: D-^Xeu, B^Xeu, and B ^ X tv. 
In what follows we discuss the implications of each of these processes in turn. 

For D —* X ev, following ref. (with minor modifications explained below), we compute 
the rate T{D ^ Xev) including order l/m?c ^^^d order Ug corrections, and compare it to 
the experimental value of BR {D ^ X e i')/td. We use 

BR{D^ ^ Xe^u) = l7.2±1.9%, r(D^) = 1.066 ± 0.023 ps. (5) 

The theoretical expression for the decay rate can be found in ref. [Q. Demanding that 
the theoretical prediction be compatible with the experimental result, we constrain the 
allowed range in the A — Ai plane. This is plotted in fig.l. The allowed region is the area 
above the curve marked with (1). A large uncertainty arises in the calculation because the 
0{a1) corrections are unknown. Effectively, these corrections set the scale at which the 
0{as) corrections should be evaluated. The curve in the figure results from the rather safe 
assumption that the scale is not above rric. If one assumes that the scale is larger than some 
value, say nic/S, then an upper bound on Ai can also be obtained. Such a bound does not 
play an essential role in our analysis, so we do not present it here. The interested reader 
may consult ref. 0. In deriving our bound we took a more conservative approach than ref. 
Pj. We take exactly the same ranges for the various input parameters (except for a^), but 
allow each of them to independently vary within their respective la bounds. To indicate the 
difference between the bounds derived in these two different ways, assuming A > 0.237 GeV 
ref. 1^ quotes Ai > — O.SGeV^ while we obtain Ai > — 0.65GeV^. Our bound is likely to 
strengthen once 0{al) corrections are calculated or if the experimental lower bound on the 
branching ratio increases. 

The decay B ^ X eu is analyzed in a similar way, namely we calculate the rate r(i? — >■ 
X eu) including order l/ml, l/ml and a^ corrections, and compare it to the experimental 
value of BR {B ^ X e v)/tb- We use for the branching ratio |^ 

Si?(5^Xei^) = 10.7 ±0.5%, (6) 



where we combined the systematic and statistical errors in quadrature, and for r^lKbp the 
most recent model independent determination from exclusive B -^ D*ii' decays | p!o[| , 



-4|-|Kfe|= 0.037 ±0.007. (7) 

1.49ps 

This procedure yields the upper bound on A corresponding to the curve marked with (2) 
in fig.l. The effect of 0{a'^) corrections is likely to strengthen our bounds (by lowering 
the effective scale for a^) but is expected to be small. The sensitivity to higher order 
corrections in the relations (^ is large. In deriving this bound we included estimates of 
these uncertainties. If the upper bound on | Vcfe| becomes lower, that would make this bound 



significantly stronger. For example, \ y^~s l^^fcl < 0.040 would yield the upper bound given 

by the solid curve in fig. 1. Similarly, an improvement in the lower bound on the branching 
ratio would strengthen our bound. 

The decay B -^ X t u has been investigated in detail in ref . [Q , where the theoretical 
expression for the branching ratio can be found, including order as and l/ml corrections. 
This decay rate, normalized to the experimentally measured rate for B ^ X eu, is a very 
reliable way of setting bounds on A and Ai. The sensitivity to higher order QCD corrections 
is minute, as the correction to the ratio of rates is much smaller than the correction to each 
of them separately 0]. In addition, the ratio is independent of the overall mf |K;fep factor. 
Taking into account the various uncertainties in the same conservative approach described 
in ref. H, and using the recent measurement |TT| with its statistical and systematic errors 



combined in quadrature 

BR{B ^ XTu) = 2.76±0m%, (8) 

we obtain a bound excluding large values of both A and Ai. This is the curve marked with 
(3) in fig.l. This bound could become much stronger if the lower bound on the B to r 
branching ratio becomes stronger. To demonstrate that, we plot (dashed line in fig.l) the 
bound that would correspond to a/2 of the above range, i.e. BR [B -^ X t u) > 2.44%. 

Each of the three experimental results that we have discussed determines an allowed hand 
in the A — Ai plane. We have not presented the upper bound from D — »• X e i/ because of 
the theoretical uncertainties discussed above. As for the lower bounds from B ^ X ev and 
B ^ X TV, aX present they are too weak to be interesting (excluding only negative values 
of A and large negative values of Ai). Therefore, we do not present them either. 

We have also investigated the implications of existing data on the inclusive rare decays 
B ^ Xs'j and B -^ Xsi^£~. These processes yield no useful constraints at present. It 
is interesting to note, however, that as the lower bound on rrit will increase and the upper 
bound on BR {B -^ Xg 7) will decrease, the resulting bounds will be almost parallel (in the 
A — Ai plane) to the bound provided hj D ^ X eu decays. For example, if experiments find 
rrit ~ 170 GeV and BR (5 — > X^ 7) < 3 x 10~^, then we will obtain the upper bound given 
by the dash-dotted curve in fig. 1. This will be important as the theoretical uncertainty in 



this process [|T2[ is much smaller than in D ^ X eu decays. 

Our results are summarized in fig.l. We learn that the two parameters A and Ai are 
significantly constrained by existing experimental data. Our most important result, an upper 
bound on A, can be easily read off from this figure: 



A<0.7GeV. (9) 

To set an upper bound on Ai, we would need to know the size of the O^a"^) corrections to 
semileptonic D decays. Assuming that the relevant value of the strong coupling constant ag 
in D decays is as{mc), ag < 2as{mc), ag < 3as{mc), we obtain the bounds Ai < 0.6 GeV^, 
Ai < 0.7 GeV^, Ai < 0.8 GeV^, respectively. The bound Ai < 0.8 GeV^ can also be obtained 
by combining the QCD sum rules prediction for A with the experimental results for BR [B — »• 
Xtu) §. 

To set a lower bound on Ai, we need a lower bound on A. We should mention a theoretical 
model-independent bound [IT3 A > 0.237 GeV (plotted with a dotted line in fig. 1). However, 



in ref. [0] it is argued that this bound is not valid. If A > 0.237 GeV, we obtain Ai > 
-0.65 GeV^ (and also Ai < 1.2 GeV^). The QCD sum rule prediction A > 0.50 GeV ^^ 
yields Ai > -0.10 GeVl 

It is instructive to compare these bounds to various model calculations. Our upper bound 
on A is hardly weaker than the allowed range of this parameter predicted by QCD sum rules, 
A = 0.57±0.07 GeV ||l|JT5[]. The QCD sum rules estimates of the parameter Ai varied strongly 



over the past few years []TB[ . These sum rules are rather unstable, thus their predictions are 



uncertain and have large errors. Recent progress in this direction |T^, however, indicates 
that Ai is likely to be negative and small in magnitude, —0.3 ^ Ai < GeV^ [|^. This is 
certainly allowed by our bounds, but then A is likely to be smaller than 0.5 GeV. A recent 
calculation of Ai in the ACCMM model resuhs Ai ~ -0.08 GeV^ H. 

These experimental constraints on A and Ai can be easily translated into bounds on the 
charm and bottom quark masses. We find that in the experimentally allowed region 

nic > 1.43 GeV , rrib > 4.66 GeV . (10) 

We do not quote upper bounds on the quark masses, both because of the uncertain status 
of the model independent lower bound on A, and because such a bound would also depend 
on assumptions about the higher order QCD corrections to semileptonic D decays. 

A precise knowledge of the bottom and charm quark masses (and of Ai) is important, 
for example, for a model independent determination of \Vcb\ from inclusive semileptonic 
B — > Xciiy decays. When phenomenological constraints, like those derived in this paper, 
or possible theoretical model-independent bounds on A and Ai will become more restrictive, 
inclusive B —>■ Xciiy decays may provide an alternative determination of | Vcb \ , comparable 
in accuracy to that from the zero recoil limit of exclusive B — > D^*Hu decays. 

The bounds for the quark masses in eq. (p!0|) have important consequences for BR {B — »• 
Xeu). The theoretical result for the semileptonic branching ratio derived by calculating 



both the semileptonic and hadronic widths, seems to constitute a puzzle at present pO 
theory predicts BR{B ^ X eu) > 0.12, in confiict with the experimental result (P). While 
we have nothing to add here to the ongoing discussion on the uncertainties in the calculation 
of the hadronic width, we would like to recall that a solution of the problem prefers low quark 
masses [^. Thus, our lower bounds on rUc and rrib point into the direction of increasing 
this problem, rather than eliminating it. 

To summarize our main results, we showed that a combination of experimental data on 
semileptonic heavy meson decays gives the upper bound A < 0.7 GeV and lower bounds on 
the heavy quark masses: m^ > 1.43 GeV and rUb > 4.66 GeV. Improved experimental results 



or theoretical calculations of 0{a1) corrections will make these constraints significantly 
stronger. 
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FIGURES 

FIG. 1. The allowed range in the A— Ai plane. The bounds are depicted by solid lines shaded on 
their excluded sides: (1) from D ^> X eu decays; (2) from B ^ X eu decays; (3) from B ^ X tv 



decays. The solid line would be the upper bound if experiment finds a/ i 49ps l^^l — 0.040. The 
dashed line would be the upper bound if experiment finds BR{B -^ Xtv) > 2.44%. The 
dash-dotted curve would be the upper bound if experiment finds BR (B ^ Xg^) < 3 x 10"*^ 
and mt ~ 170 GeV. 
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